CL_rho <-
function(normH, gamma){
	gamma * normH * min(1, normH)
}

